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ABSTRACT 

o 

' Context. Since the discovery of the first transiting extrasolar planet, transit timing has been recognized as a powerful method to 

discover and characterize additional planets in these systems. However, the gravitational influence of additional planets is not the only 
expected source of transit timing variations. 

Aims. In this work, we derive the expected detection frequency of stellar companions of hot-jupiter transiting planets host-stars, 
detectable by means of transit timing analysis. Since roughly half of the stars in the solar neighborhood belong to binary or multiple 
stellar systems, the same fraction of binary systems may be expected to be present among transiting planet-host stars, unless planet 
formation is significantly influenced by the presence of a stellar companion. Transit searches are less affected by the selection biases 
against long-period binaries that plague radial velocity surveys. 

Methods. We considered the frequency and the period, mass ratio and eccentricity distributions of known binary systems in the solar 
neighborhood, and estimated the fraction of transiting planet-hosts expected to show detectable transit timing variations due to the 
P_] , light travel time effect in a binary stellar system, in function of the time since the discovery of the planet. 

Results. If the frequency of binaries among hot-jupiter planet host stars is the same as determined in the solar neighborhood, after 5 
years since the discovery of a sample of transiting planets 1.0% ± 0.2% of them have a probability > 99% to present transit timing 
variations > 50 sec induced by stellar binarity, and 2.8% ± 0.3% after 10 years, if the planetary and binary orbits are coplanar. 
Considering the case of random inclinations the probabilities are 0.6% ± 0.1% and 1.7% ± 0.2% after 5 and 10 years respectively. Our 
estimates can be considered conservative lower limits, since we have taken into account only binaries with periods P > 5 • 10 3 days 
(a > 6 AU). Our simulations indicate that transit timing variations due to the light travel time effect will allow us to discover stellar 
companions up to maximum separations equal to a ~ 36 AU after 5 years since the discovery of the planet {a ~ 75 AU after 10 years). 
Conclusions. Comparing the results of the observations with the above predictions allows to understand if stellar companions at 
critical separations (< 100 AU) are favoring or hindering the formation of hot-jupiter planets. Comparing the results of transit timing 
detections with those obtained by other complementary methods results in a more complete determination of stellar multiplicity 
around transiting planet-host stars. Moreover, transit timing analysis allows us to probe stellar multiplicity at critical separations 
(< 100 AU) around stars located in different regions of the Galaxy, and not just in the solar neighborhood. 
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c^*) 1 . Introduction planets is unlikely in equal mass binaries with semi-major axis 

f^) \ ~ 50 AU, both by means of the disk instability (Boss 1997) and 

t-H , The frequency of binary stars and multiple stellar systems the core accretion mechanisms (Pollack et al. 1996). On the con- 

. around solar-like stars in the solar neighborhood has been ex- [mry Boss {2 006), showed that in the context of the disk insta- 

. £ ; tensively studied in the past. Duquennoy & Mayor (1991), from b0ity scen ario, the presence of a close-by stellar companion may 

V , a sample of 164 primary G-dwarf stars analyzed during almost in fact ^gger dumps formation leading to giant planets. Along 

u ■ 13 yr with the CORAVEL spectrograph, obtained that the ratios the same lines Duchene (2010) suggested that in tight binaries 

Cd ; of single:double:triple:quadruple systems are 51:40:7:2 respec- (fl < m AU) massive planets ^ formed by disk instabi i ity at 

tively. This fraction of binaries should be present also in a sam- the same rate as less massive planets in wider binaries and single 

pie of planet host stars if planets do form in any kind of binary stars 

systems and there are no selection biases applied to define the According to Marzari et al. (2005), independently from the 

sample of stars where the planets are searched for. planet formation mec hanism, tidal perturbations of the com- 

From a theoretical point of view the problem of the forma- panion star i nflue nce both the onset of instability and the fol- 

tion of giant planets in close binary systems is still largely de- lowing chaotic evolution of protoplanetary disks. In particu- 

bated. Truncation and heating of circumstellar protoplanetary lar several studies have po i nt ed out that the gravitational influ- 

disks are expected to occur in close binary systems (Artymowics ence of stellar compan i ons on the dynamics of planetary sys- 

& Lubow 1994; Nelson 2000). However there is no general con- t ems becomes significant at separations < 100 AU (e.g. Pfhal & 

sensus on which consequences these processes have on planet Muterspaugh 2006; Desidera & Barbieri 2007; Duchene 2010). 

formation. Nelson (2000) showed that the formation of giant DeS pite that, four planetary systems have been discovered in bi- 

naries with separations around 20 AU: Gamma Cephei (Hatzes 

Send offprint requests to: M. Montalto et al. 2003), Gliese 86 (Queloz et al. 2000; Els et al. 2001), HD 
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41004 (Santos et al. 2002), and HD196885 (Correia et al. 2008). 
These observational results are clearly challenging our current 
knowledge of planet formation and evolution in binary systems. 

Several imaging surveys have successfully identified stel- 
lar companions to planet-host stars (e.g. Patience et al. 2002; 
Luhman & Jayawardhana 2002; Chauvin et al. 2006; Mugrauer 
et al. 2007; Eggenberger et al. 2007; Eggenberger & Udry 
2007). At present, giant planets around wide binary systems 
appear as frequent as planets around single stars (Bonavita & 
Desidera 2007), suggesting that wide binaries are not signif- 
icantly altering planet formation processes. However, there is 
a marginal statistical evidence that binaries with separations 
smaller than 100 AU may have a lower frequency of planetary 
systems than around single stars. Eggenberger et al. (2008) esti- 
mated a difference in the binary frequency between carefully se- 
lected samples of non planet-host and planet-host stars ranging 
between 8.2% + 5.0% and 12.5% ± 5.9%, considering binaries 
with semi-major axis between 35 AU and 250 AU, further point- 
ing out that this difference seems mostly evident for binaries 
with a < 100 AU. These results are obtained using samples of 
stars hosting radial velocity discovered planets, for the obvious 
reason that Doppler spectroscopy has been so far the most suc- 
cessful planet detection method, providing then the largest sam- 
ple of planets from which statistical conclusions can be drawn. 
Nevertheless, planets discovered with this technique are known 
to be adversely selected against close binaries. Doppler spec- 
troscopy is complicated by light contamination of stellar com- 
panions, and blending of spectral lines. This implies that radial 
velocity surveys typically exclude known moderately and close 
binaries from their target lists. As a consequence the occurrence 
of planetary systems in binaries, and in particular at critical sep- 
arations (< 100 AU), is still poorly constrained by the obser- 
vations. While the use of a control sample of stars, proposed 
by Eggenberger et al. (2007, 2008), is expected to mitigate the 
impact on the results given by the bias against close binaries 
applied by radial velocity surveys, other samples of planet-host 
stars and different techniques may be useful to probe the fre- 
quency of planets in binary stellar systems. 

While radial velocity planet searches tend to exclude binaries 
from their target lists, transit searches do not apply a priori se- 
lection criteria against binaries. Transiting planets are routinely 
searched among all stars, and then subjected to follow-up pho- 
tometric and spectroscopic analysis. Photometric analysis aims 
mainly at ruling out grazing eclipsing binary stellar systems 
which manifest themselves by means of markedly V-shaped 
eclipses, the presence of secondary eclipses, color changes dur- 
ing the eclipses and light modulations with the same periodicity 
of the transiting object. Spectroscopic analysis is then used to 
further rule-out giant stars primaries and the more complicated 
scenarios involving hierarchical triple systems with an eclips- 
ing binary stellar system, and blends with background eclipsing 
binaries (Brown 2003). However, these follow-up analysis do 
not eliminate planets in binary stellar systems. Several transiting 
planets are already known members of binaries (e. g. Daemgen 
et al. 2009), and others have suspected close companions as in- 
dicated by the presence of radial velocity and transit timing vari- 
ations (Winn et al. 2010; Maxted et al. 2010; Queloz et al. 2010; 
Rabus et al. 2009). Moreover for transiting planets we have 
a firm constraint on the planetary orbital inclination (which is 
close to 90°). The Rossiter-McLaughlin effect (Rossiter 1924; 
McLaughlin 1924) can be used to probe the sky-projected an- 
gle j6 between the stellar rotation axis and the planet's orbital 
axis. By transforming the projected angle /3 into the the real spin- 
orbit angle if/ using a statistical approach and the entire sample 



of planets with Rossiter-McLaughlin measurements, Triaud et 
al. (2010) derived that most transiting planets have misaligned 
orbits (80% with tff > 22°), and that the histogram of pro- 
jected obliquities closely reproduces the theoretical distribution 
of iff using Kozai cycles and tidal friction (Fabrycky & Tremaine 
2007). Since type I and II migration are not able to explain the 
present observations, the indication is that the Kozai mechanism 
is the major responsible of the formation of hot-jupiter planets. 
In this case, we should expect that most hot-jupiter planets have 
stellar companions. The discovery of close stellar companions to 
transiting planets systems is then of primary importance, since it 
would constitute a strong proof in favor of the Kozay cycles and 
tidal friction mechanism. 

Transiting planets host stars then constitute an interesting 
sample of objects where to look for additional distant compan- 
ions. The presence of a stellar companion around these stars 
can be inferred at least by means of four independent and com- 
plementary techniques: transit timing variations, radial velocity 
drifts, direct imaging and IR excess. 

Deriving the expected frequency of transiting planet host 
stars in binary stellar systems detectable by each one of the 
above mentioned techniques is then important, since compar- 
ing the results of the observations with the predictions we can 
better understand which influence close binary systems have on 
planet formation and evolution. If, for example, the existence of 
hot-jupiters is connected to the presence of close-by stellar com- 
panions, we should expect to derive a higher binary frequency 
around stars hosting these planets, with respect to the frequency 
of binaries observed in the solar neighborhood. If, on the con- 
trary, the presence of close-by stellar companions strongly pre- 
vents the existence of planets, we should expect a lower fre- 
quency. In this paper we focus our attention on transit timing 
variations (TTVs) induced by stellar binarity. Future contributes 
will account for the other techniques. 

In particular here we derive the expected frequency of tran- 
siting planets in binary systems detectable by TTVs. We define 
the detection frequency {f d e i) as the fraction of transiting plan- 
ets expected to show detectable TTVs induced by stellar binarity 
over some fixed timescales, when the only source of TTVs un- 
der consideration is the light travel time effect in binary systems. 
The presence of an additional stellar companion around a tran- 
siting planet-host star, should induce TTVs even if we neglect 
perturbing effects, because of the variable distance of the host 
star with respect to the observer in the course of its orbital rev- 
olution around the barycenter of the binary stellar system. This 
motion induces TTVs affecting the observed period of the tran- 
siting object, and consequently the ephemerides of the transits 
(e.g. Irwin 1959). 

Transit timing allows detenction of close stellar compan- 
ions around more distant planet-host stars than direct imag- 
ing. Accurate transit timing measurements are achievable also 
for planets discovered around faint and distant planet-hosts, by 
means of a careful choice of the telescope and the detector (e. g. 
Adams et al. 2010). On the contrary, the distance of the planet- 
host constitutes a limit for direct imaging detection of stellar 
companions. Using VLT/NACO, and targeting solar type close- 
by stars (~ 10 pc), a companion with a mass of 0.08 M (then 
just at the limit of the brown dwarf regime) can be detected at the 
3-cr limit at a projected separation of 0.3 arcsec (e.g. Schnupp et 
al. 2010; Eggenberger et al. 2007), which corresponds to 3 AU. 
If, however, the star is located at a distance > 333 pc, direct 
imaging can probe only separations > 100 AU. Then transit tim- 
ing allows us to probe stellar multiplicity at critical separations 
< 100 AU (as demonstrated in this work) around more distant 
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samples of transiting planet host-stars with respect to what can 
be done by direct imaging, giving the opportunity to probe stel- 
lar multiplicity around targets located in different regions of the 
Galaxy. Moreover, while direct imaging is more efficient in de- 
tecting companions in face-on orbits, transit timing (and Doppler 
spectroscopy) is more efficient in the case of edge-on systems. 

This paper is organized as follows: in Sect. [2] we review the 
known properties of binary stellar systems in the solar neigh- 
borhood; in Sect. [3j we discuss the light travel time effect of 
transiting planets in binaries; in Sect [4] we describe the Monte 
Carlo simulations we did to constrain the frequency of transiting 
planet-host stars presenting detectable transit timing variations 
induced by binarity over some fixed timescales; in Sect. [5] we 
discuss the results of our analysis; in Sect. [6j we summarize and 
conclude. 



2. Properties of multiple stellar systems 

From their study of multiple stellar systems in the solar neigh- 
borhood, Duquennoy & Mayor (1991) derived the following 
properties for binary stellar systems with mass ratios q > 0.1: 
1) the orbital period distribution can be approximated by: 



fn m r -dog P- log P) 2 
f(logP) = Cexp — 



2cr 



(1) 



logP 



where log P = 4.8, cri og p — 2.3, and P is in days; (2) binaries 
with periods P > 1000 days (which is the range of periods in 
which we are interested, see below) have an observed eccentric- 
ity distribution which tends smoothly toward h(e) = 2e; (3) the 
mass-ratio (q = ni2/mi) distribution can be approximated by: 



g(q) = kexp 



-(q-q) 2 
2o% 



(2) 



where q = 0.23, cr q = 0.42, k = 18 for their G-dwarf sample. 

In the following we will consider only binaries with periods 
P > 5 ■ 10 3 days. This limit is due to the need to minimize 
perturbing effects, which have been explicitely neglected in our 
analysis, as explained in Sect. [3] 

Using the period distribution of Duquennoy & Mayor 
(1991), we have that ~68% of all binary systems are expected 
to have P > 5 ■ 10 3 days and, considering that the frequency 
of binary stellar systems in the solar neighborhood is 40% (as 
reported in Sect. Q] excluding multiple stellar systems), we have 
that the expected frequency of binaries in our period range is 
equal to 68% ■ 40% - 27%. 

The aim of the rest of this paper is to establish how many of 
these systems should appear as detectable transit timing sources 
over a timescale of at most 10 years since the discovery of the 
transiting planets around the primary stars of these systems. 



3. The light travel time effect 

In the following we consider the case of a planet orbiting the pri- 
mary star of a binary stellar system. This configuration is usually 
called S-type orbit. We will also assume that the planet is re- 
volving much closer to its parent star than the companion star 
(Pbin » Ppi where f/„„ is the binary orbital period, and P p \ is 
the planet orbital period). This is a reasonable assumption given 
the range of binary periods we are considering (Pbin > 5 ■ 10 3 
days, see above), and that most transiting planets are typically 



hot-jupiters with a period of only a few days. Under this as- 
sumption, and given also that our adopted observing window 
(10 years) is smaller than the shortest period we considered, the 
perturbing effects of the secondary star on the orbit of the inner 
planet can be neglectedQ The three body problem can be splitted 
in two independent two body problems: the motion of the planet 
around the host star (which can be reasonably assumed coinci- 
dent with the barycenter of the planetary system given that the 
mass of the planet is much smaller than the mass of the star), 
and the motion of the host star around the barycenter of the bi- 
nary system. Even neglecting perturbing effects, we expect tran- 
sit timing variations to be present, due to the light travel time 
effect, as described in Sect. Q] The transiting planet can be used 
as a precise clock unveiling the presence of the additional star in 
the system. 

We consider a reference system with the origin in the 
barycenter of the binary system. The Z axis is aligned along 
the line of sight, the X axis along the nodal line defined by the 
intersection between the binary orbital plane and the plane of 
the sky, and the Y axis consequently assuming a right-handed 
Cartesian coordinate system. The XT-plane is tangent to the ce- 
lestial sphere. 

The distance (z) between the host star and the barycenter of 
the binary system, projected along the observer line of sight is 
given by (e. g. Kopal 1959): 



z = 



fli(l -e 2 ) 
1 + ecos(f) 



sin(u) + /) sin(i), 



(3) 



where the orbital elements are relative to the barycentric orbit 
of the host star. Using the equation of the center of mass a\ = 
ni2al(m\+ mi), the Kepler law and dividing by the speed of light 
(c), gives: 



0(0 = 



( G )l/3 p2/3 (1 _ g 2) m sin ^ sin ^ + y W ) 



c (27r) 2 / 3 (m! + m 2 ) 2/3 1 + e cos(f[t]) ' 



(4) 



where / is the true anomaly, e is the eccentricity of the orbit, and 
i is the inclination with respect to the plane of the sky, m\ and 
m.2 are the masses of the primary and of the secondary, co is the 
argument of the pericenter, P is the period, G the gravitational 
constant, t the epoch of observation, and we just recall that the 
orbital elements are relative to the binary orbit. 

Eq. 4 gives the time necessary to cover the distance projected 
along the line of sight between the host star and the barycenter 
of the binary system at the speed of light. We can also say that 
Eq. 4 gives the difference between the observed ephemerides 
of the transiting planet once including the light time effect and 
the ephemerides obtained considering the intrinsic period of the 
planet . 

The rate of change over time of the observed planet period 
with respect to the intrinsic planet period is given by the deriva- 
tive of Eq. 4 over time: 



e(t) = 



(2tt) 1/3 (G) 1/3 to 2 sin(i) 



cP 1 ' 3 VT^Cmi +m 2 ) 2 / 3 



{cos(u+f[t])+e cos((jS) j.(5) 



The procedure usually adopted by observers to determine the 
period of a transiting object is based on the determination of the 



1 In particular secular perturbations, see e.g. Kopal (1959) for a dis- 
cussion of perturbing effects of a third body on a close eclipsing binary. 

2 We are assuming that the barycentric radial velocity of the binary 
stellar system is y = 0. 
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time interval between two or more measured transits at an epoch 
to. Then e(fo) = £q can be seen as the difference between the ob- 
served and the intrinsic period at the observed time to due to light 
time effect; to avoid accumulation of this difference over time in 
the ephemeris it needs to be subtracted. Fig.Q](upper figure, up- 
per panel, dotted line) shows this accumulation with time. The 
difference between the observed and the predicted ephemerides 
(O - C) of successive transits at another generic epoch (charac- 
terized by the binary true anomaly / and the epoch t — t q + A f), 
is then given by: 



[0-C](t) = 0(t) - O(t ) - e (t - t ). 



(6) 



Imposing that the O - C residual is larger than a given detection 
threshold r is equivalent to solve the following disequation: 



\O(t + At) - O(to) - e (At) \ - r > 0. 



(7) 



In Fig. Q] (upper figure), we present an illustrative example. We 
consider a binary system where the relevant orbital elements 
were fixed at the median values derived by Duquennoy & Mayor 
(1991): P bin =172.865 years, e = 0, q = m 2 /m 1 = 0.23, to = 0, 
i = 90°. We assume that a planet is orbiting the primary star of 
the binary system, and that the period of the planet was mea- 
sured at the epoch to (assumed coincident with the origin of the 
x-axis in Fig.[TJ. After five years the transit timing variation is 46 
sec, and after ten years it is 178 sec, as shown in the inner plot 
of Fig. Q] The binary we have considered is one of the typical 
binaries in the solar surrounding. We also observe that the O — C 
diagram is in general not periodic as evident already from Fig.Q] 

4. Simulations 

In this Section we describe the orbital simulations we did to con- 
strain the expected frequency of transiting planet-host stars that 
should present a detectable transit timing variation induced by 
stellar binarity over a timescale of 5 years and 10 years since the 
discovery of the planet around the primary star of the system. 
The detectability threshold was fixed considering the results ob- 
tained by transit timing searchers with ground based telescopes. 
Rabus et al. (2009), observing the planetary system TrES-1 with 
the IAC80 cm telescope obtained mean precisions of 18.5 sec, 
which is also in agreement with the theoretical equation given 
by Doyle & Deeg (2004). On the basis of that, and using sev- 
eral other literature results, they were already able to claim the 
presence of a linear trend equivalent to 48 sec (obtained from 
their best-fit parameter) over the period of 4.1 yr spanned by 
the entire sample of the observations. Winn et al. (2009) mea- 
sured two transits of the giant planet WASP-4b with the Baade 
6.5 m telescope obtaining mid-transit times precisions of ~6 sec. 
Typically it can be assumed that small or moderate groundbased 
telescopes can reach precisions < 20 sec, while large telescope 
may reach better than 10 sec precision. Given these results we 
considered that a TTV detection threshold equal to r =50 sec 
can be reasonably applied in our analysis. 

We randomly chose the period, mass ratio, and eccentricity 
of the binary considering the probability distributions presented 
in Sect. [2] The mass of the primary star was fixed at 1 M , since 
we are considering binarity among typical solar type stars. The 
argument of the pericenter was instead randomly chosen using a 
uniform probability distribution. The inclination was either fixed 
to 90° (to consider the case of coplanar orbits) or randomly cho- 
sen using a uniform probability distribution. For each orbit (char- 
acterized by P, e, q, u>, i) we numerically solved the disequation 
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Fig. 1. Upper figure, upper panel: the black solid curve repre- 
sents the difference between the observed transit ephemerides 
of the planet once including the light time effect and the 
ephemerides obtained once considering the intrinsic period of 
the planet (Eq. 4, see text). The dotted line represents the accu- 
mulation over time of the difference between the observed pe- 
riod of the planet and its intrinsic period if the planet period is 
measured at the epoch to coincident with the origin of the x-axis, 
and it is then assumed constant. Upper figure, lower panel: the 
O - C diagram (Eq. 6). The inner plot shows a close-up view 
of the O - C diagram during the first 10 years after the period 
determination. Units are the same of the large plot. We assumed 
P Mn =172.865 years, e = 0, q = m%jm\ = 0.23, tj = 0, i = 90°. 
Lower figure: Graphical representation of disequation 7 in func- 
tion of to, the epoch of the planetary period determination, for 
the case of the orbit considered in the upper figure. The horizon- 
tal black solid lines denote the time intervals where the condition 
for transit timing detectability is met. The value of r is 50 sec, 
and the timescale A T is 10 years (see text for details). 
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Fig. 2. Distributions of the detection probabilities (Pder) of tran- 
sit timing variations of transiting planets in binary stellar sys- 
tems, relative to our sample of 10 6 simulated binary orbits. The 
simulations are performed in function of different time intervals 
since the discovery of the planets (At = 5 years or 10 years), 
and in function of different assumptions on the inclination of the 
binary stellar system orbit, as indicated by the labels in the dif- 
ferent panels. 



(7) in function of ?o (the epoch of the planetary period deter- 
mination). We subdivided the period P in 10000 equal intervals 
of time and evaluated disequation (7) at the extremes these in- 
tervals. Then we isolated the intervals in which disequation (7) 
changed sign, and using the secant method imposing a threshold 
for the convergence equal to 0. 1 sec we obtained the roots of the 
correspondent equation. Then we determined the intervals A to 
where disequation (7) was satisfied. Summing up together these 
intervals of time and dividing by the period of the binary gave 
the probability to observe the requested transit timing variation 
for that fixed orbit over the given timescale (At - t—to either 5 yr 
or 10 yr in our simulations) assuming to determine the period of 
the transiting planet in correspondence of a random orbital phase 
of the binary. In such a way we assigned to each simulated orbit 
a transit timing detection probability (Pdei)- In Fig. Q~|(lower fig- 
ure), we show the graphical representation of disequation (7) in 
function of fo, the epoch of the planetary period determination, 
for the case of the orbit considered in Fig. Q] (upper figure), as- 
suming a transit timing threshold t — 50 sec, and a timescale 
of 10 years. The horizontal black solid lines denote the time 
intervals where disequation (7) is satisfied. We performed 100 
runs of 10000 simulations each, calculating the mean detection 
probabilities and their 1-cr uncertainties, as reported in the next 
Section. 



5. Results 

The final detection probability histograms are shown in Fig. [2] 
obtained from the entire sample of 10 6 simulated orbits. While 
most orbits imply null detection probabilities over the assumed 



timescales, in each one of the different situations we considered, 
the histograms present a probability tail extended toward large 
detection probabilities. The number of orbits having Pdei > 99% 
is equal to 3.9% ± 0.2% after 5 years since the period determina- 
tion (10.4% ±0.3% after 10 years) for the case of coplanar orbits, 
and is equal to 2.2% ±0.1% after 5 years (6.2% ± 0.2% after 10 
years) for the case of random inclinations. The fact that the his- 
tograms are extended toward large probabilities is a consequence 
of the period distribution and the adopted timescales. Orbits hav- 
ing Pdei > 99% have also periods smaller than P < 1.67 ■ 10 5 
(considering a timescale of 10 years), which means that transit 
timing can allow discovery of stellar companions up to separa- 
tions equal to a ~ 75 AU after 10 years since the discovery of 
the planet (a ~ 36 AU after 5 years). 

Then, considering the observed frequency of binaries in the 
solar surrounding with periods P > 5 ■ 10 3 days (27%) and the 
case of coplanar orbits, after 5 years since the discovery of a 
sample of transiting planets 1.0% ± 0.2% transiting planet host- 
stars will have a probability Pdet > 99% to present detectable 
(> 50 sec) transit timing variations induced by stellar binarity, 
and 2.8% ± 0.3% after 10 years. Considering the case of ran- 
dom inclinations the expected frequencies (fdet) are 0.6% ±0.1% 
and 1.7% ± 0.2% after 5 and 10 years respectively. These results 
are summarized in Table Q] Our estimates can be considered a 
conservative lower limit, since we have excluded binaries with 
periods P < 5 ■ 10 3 days. 

6. Conclusions 

In this paper we have investigated how known transiting extra- 
solar planets can be used to constrain the frequency of multiple 
stellar systems among planet-host stars. The presence of a stellar 
companion in these systems is expected to induce transit timing 
variations of the transiting planets even once perturbing effects 
are neglected, due to the orbital revolution of the primary around 
the barycenter of the binary stellar system. 

If the frequency of binaries among planet-host stars is the 
same as determined in the solar neighborhood, after 5 years since 
the discovery of a sample of transiting planets 1.0% ± 0.2% of 
them have a probability > 99% to present a transit timing varia- 
tions > 50 sec induced by stellar binarity, and 2.8 ±0.3% after 10 
years if the planetary and binary orbits are coplanar. Considering 
the case of random inclinations the probabilities are 0.6% ±0.1% 
and 1.7% ± 0.2% after 5 and 10 years respectively. Our results 
have been obtained assuming a binary period P > 5 ■ 10 3 (a > 6 
AU). Moreover, we derived that we can expect to discover stel- 
lar companions of transiting planets host stars up to a maximum 
separations a ~ 75 AU after 10 years since the discovery of a 
planet (a ~ 36 AU after 5 years). 

A final comment is necessary to mention that TTVs may 
have several different origins among which perturbing effects 
caused by additional planets or moons, secular precession due 
to general relativity, stellar proper motion, the Appelgate ef- 
fect (Agol et al. 2005; Miralda-Escude 2002; Nesvorny 2009; 
Heyl & Gladman 2007; Ford & Holman 2007; Simon 2007; 
Kipping 2009a; Kipping 2009b; Pal & Kocsis 2008; Rafikov 
2009; Watson 2010), and binarity of the host is one of them. 
However, transit timing variations induced by binarity are ex- 
pected to produce long-term trends, and in particular they should 
be associated also with radial velocity drifts of the host star. 
Transit timing searchers should also follow-up spectroscopically 
their targets, since a transit timing variation associated with a ra- 
dial velocity variation will be very likely the signature of bina- 
rity. 
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Table 1. Expected frequency (fdet) of binary stellar systems de- 
tectable by means of transit timing variations in function of the 
timescale (A T) since the discovery of the planet, and of different 
assumptions on the inclination of the binary stellar systems. 



i = 90° 


AT 


fdet 




(yr) 


(%) 




5 


1.0 ±0.2 




10 


2.8 ±0.3 
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http://www.superwasp.org/documents/triaud2010_rossiter.pdf 
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0°<;<90° AT f de , 

(yr) (%) 

5 0.6 ±0.1 

10 1.7 ±0.2 
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